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An optical realization of the tunneling dynamics of two interacting bosons in a double-well po- 
tential, based on light transport in a four-core microstructured fiber, is proposed. The optical 
setting enables to visualize in a purely classical system the entire crossover from Rabi oscillations 
to correlated pair tunneling and to tunneling of a fragmented pair in the fermionization limit. 
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I. INTRODUCTION. 

Light transport in engineered optical waveguides has 
provided a fascinating and experimentally accessible 
framework to visualize in a classical setting many univer- 
sal coherent quantum phenomena generally encountered 
in condensed-matter or matter- wave systems [TJ [5] . This 
has lead to the prediction and observation of a wide va- 
riety of classic optics analogues of single-particle nonrel- 
ativistic and even relativistic phenomena, such as Bloch 
oscillations and Zener tunneling [TJ [3] , dynamic localiza- 
tion [3] , Anderson localization [5] , coherent destruction of 
tunneling j2], Zeno dynamics [7], adiabatic stabilization 
[8], and Zitterbewegung [9]. Since photons do not inter- 
act, it is a common belief that, as opposed to other quan- 
tum systems such as cold atoms or trapped ions (see, e.g., 

[10) ). the use of photonics as a model system for quan- 
tum physics carries the intrinsic drawback of being lim- 
ited to visualize single-particle phenomena, missing the 
possibility to simulate the richer physics of interacting 
many-particle quantum systems. A paradigmatic exam- 
ple of many-body physics is found in quantum tunneling 
of bosons in a double well potential, the so-called bosonic 
junction |llH17j . For a relatively large number and 
weakly interacting bosons, this has led to the observa- 
tion of Josephson oscillations and nonlinear self-trapping 
of bosons above a critical interaction strength, as de- 
scribed by a standard Bosc-Hubbard model or by cou- 
pled mean-field equations in the Gross-Pitaevskii limit 

[11) . A simple optical realization of the bosonic junction 
in such a limiting case is based on light tunneling be- 
tween two coupled nonlinear waveguides |18| . However, a 
richer dynamical scenario has been recently predicted to 
occur for tunneling of few and strongly correlated bosons 
[TFjI [TO] , covering the full crossover from weak interac- 
tions to the fermionization limit of the Tonks- Girardeau 
gas [20] . In particular, the tunneling dynamics of two 
bosons in a one-dimensional double well shows a transi- 
tion from Rabi oscillations, in the absence of interaction, 
to correlated pair tunneling and further to fragmented- 
pair tunneling as the interaction strength is increased 
[T5] . As few-body counterparts of the self-trapping tran- 
sition and correlated pair tunneling in a bosonic junc- 
tion have been reported in recent experiments [141 12 1] , 
an observation of the rich two-boson tunneling dynam- 
ics up to the fermionization limit |15j is still missing. In 



this article it is shown that such a two-boson tunneling 
dynamics can be realized in a classical optical setting 
based on four-core guiding dielectric structure, in which 
the electric field propagation along the guide mimics the 
quantum mechanical evolution of the two-particle wave 
function. The paper is organized as follows. In Sec. II, 
the quantum-optical analogy between light propagation 
in a four-core microstructured fiber and the dynamics of 
two interacting bosons in a double well is outlined. In 
Sec. Ill, a detailed analysis of the tunneling dynamics is 
presented, and the entire crossover from Rabi oscillations 
to correlated pair tunneling and to tunneling of a frag- 
mented pair in the fermionization limit is explained on 
the basis of the coupling among the various modes sus- 
tained by the fiber cores. Finally, in Sec. IV the main 
conclusions are outlined. 



II. QUANTUM-OPTICAL ANALOGY 

Let us consider a weakly-guiding dielectric structure 
with a refractive index n(xi,X2), which varies in the 
transverse {x\,X2) plane but remains invariant along the 
axial direction z. In the paraxial approximation, prop- 
agation of monochromatic light waves is described by a 
Schddinger-type wave equation for the electric field en- 
velope tp [5] 

where A = \/(2ir) is the reduced wavelength of photons, 
V(xi,X2) — n s — n(x\,X2) is the optical potential, and 
n s is the substrate refractive index. The normalization 
condition /_ dx \dx2\ip \ 2 — 1 will be assumed in the 
following. Previous quantum-optical analogies have gen- 
erally viewed the paraxial wave equation (1) as formally 
equivalent to the Schrodinger equation for a single parti- 
cle of mass n s in a two-dimensional potential V(x\,X2), 
in which the temporal evolution of the quantum parti- 
cle is mapped into the spatial light evolution along the 
axial direction z and the Planck's constant is replaced 
by the reduced wavelength of photons (see, for instance, 
[5]). However, whenever the potential V has the form 

V(xx,x 2 ) = V w (x 1 ) + V w (x 2 ) + V mt (\x 1 -x 2 \), (2) 
where V w (x) is an arbitrary one-dimensional potential 
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FIG. 1. (color online) (a) Profiles of the double well poten- 
tial V w (x) (solid curve) and of the repulsive potential Vmt{x) 
(dashed curve) for parameter values a — 4.5 /im, w = 3 /jm, 
D x — 1 /j,m, Wi = 0.5 /im, D xi = 0.2 /im, Ani = 0.003, 
and An2 = 0.002. The corresponding two-dimensional opti- 
cal potential V(xi,X2) is shown in (b). (c) Schematic of the 
guided modes supported by the four core regions involved in 
the tunneling dynamics and their couplings. 



and Vi n t(x) is a short-range potential, Eq.(l) can be re- 
garded as the optical analogue of the Schrodinger equa- 
tion for two particles with the same mass n s in a one- 
dimensional potential V w , which interacts via the poten- 
tial Vi n t- If the optical structure is excited at the z = 
input plane by a beam satisfying the symmetry constraint 
ip(xi, X2, 0) = ip(x2,Xt,0), the wave function ip remains 
symmetric along the propagation, and Eq.(l) thus de- 
scribes the evolution of two interacting identical bosons. 
Therefore, if we assume for V w a double well shape and for 
Vi n t a short-range repulsive potential, our optical system 
realizes a classic wave optics analogue of the two-boson 
junction recently studied in Ref . [T5] . In our optical sys- 
tem, we assume for V w (x) a double well of the form [6] 
V w = —Ani[g(x — a) + g(x + a)], where g(x) = [erf ((a; + 
w)/D x ) - erf((x - w)/D x )]/{2erf(w/D x )} is the well 
shape, 2a is the distance between the two wells, Arc-i > 
is the peak index change that defines the well depth, and 
2w is the well width. For the repulsive potential, we 
assume a similar functional form Vint — Ari2[erf((a: + 
Wi)/D xi ) - erf((x - Wi)/D xi )]/[2eri {w t /D xi )], with w, 
and D x i much smaller than w and D x , respectively. 
The refractive index change An 2 > measures the 
strength of the interaction, Ari2 = corresponding to 
non- interacting bosons. Typical shapes of V w {x), Vi n t(x) 
and of the resulting two-dimensional potential V(xi, x%) 
[Eq.(2)] are shown in Figs. 1(a) and (b). Note that the 
resulting potential V in the (xi,x 2 ) plane defines four 
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FIG. 2. (color online) Numerically-computed behavior of (a) 
the percentage of bosons in the right well pr, and (b) of the 
pair probability P2 versus propagation distance, for increasing 
values of the particle interaction strength, measured by the in- 
dex change Ari2. Curve 1: Ari2 = (non-interacting bosons); 
curve 2: An2 = 0.5 x 10~ 3 ; curve 3: Ari2 = 1.5 x 10~ 3 ; curve 
4: Ari2 = 15 x 10~ 3 . Curves 2 and 3 correspond to the corre- 
lated pair tunneling regime, whereas curve 4 corresponds to 
tunneling of a fragmented pair. 

higher-index guiding regions, i.e. four waveguides, de- 
noted by I-IV in Fig. 1(b), which are evanescently cou- 
pled. Such a four-core guide could be realized, for exam- 
ple, with the technology of microstructured fibers [22] , in 
which a preform with the desired geometrical and refrac- 
tive index features is first manufactured. For example, 
using a cladding region made of fused silica, the struc- 
ture of Fig. 1(b) could be realized by assembling different 
regions of fused silica with different GeC»2 doping con- 
centrations. 



III. TUNNELING DYNAMICS 

The main features of the tunneling dynamics of two 
bosons in a double-well potential are captured by ana- 
lyzing the evolution of the percentage of bosons in the 
right well pr{z) and the the pair (or same-site) boson 
probability P2(z), which are defined by |15] 

/•oo /*oo 

Pr(z) = I dxx dx 2 \ip\ 2 (3) 

JO J -00 

p 2 0)= / dx 1 dx 2 \ip\ 2 + / dx 1 dx 2 \M 2 ( 4 ) 

In our optical setting, pr{z) and P"2,(z) simply correspond 
to the fractional light power trapped in waveguides I and 
IV, and in waveguides I and III, respectively. A typical 
evolution of Pr(z) and p 2 {z), as obtained by numerical 
integration of Eq.(l) for increasing values of the inter- 
action strength An.2, is shown in Fig. 2. Parameter val- 
ues used in simulations are A = 633 nm and n s — 1.45. 
In each simulation, the structure is excited at z = in 
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the fundamental mode of the guide I, which corresponds 
to have initially the two bosons in the right-side well in 
the lowest energy state. The scenario shown in Fig. 2 re- 
produces the transition from uncorrelated tunneling to 
pair tunneling and fragmented tunneling in the fermion- 
ization limit, predicted in Ref.|15j. For non-interacting 
bosons (curve 1), the atoms simply Rabi oscillate back 
and forth between both wells, and they tunnel indepen- 
dently. As a small correlation is introduced (curve 2), 
both atoms tend to remain in the same well in the course 
of tunneling, i.e they tunnel as pairs. Such a dynamical 
behavior, which was observed in |16| and referred to as 
second-order tunneling, can be simply explained in the 
framework of a standard two-site Bose-Hubbard model, 
the optical simulation of which was recently proposed in 
the Fock space using waveguide arrays |23| . However, the 
standard Bose-Hubbard model fails to predict the tun- 
neling regimes at strong interaction and the transition to 
the fermionization limit. Indeed, at a larger interaction 
(curve 3), tunneling tends to be inhibited, which is the 
few-body signature of the self-trapping phenomenon of 
many bosons in the mean-field limit. Remarkably, at 
stronger interaction and near the fermionization limit 
(curve 4), tunneling is again allowed, and a fast oscil- 
lation of pr(z) is superimposed to the slower tunneling 
cycle. This basically corresponds to fragmented-pair tun- 
neling at the Rabi frequency predicted in Ref. [15] . Cor- 
respondingly, p 2 (z) passes through just about any value 
from 1 (fragmented pair) to small values (near complete 
isolation). A detailed explanation of such a rich tunnel- 
ing scenario requires an inspection of the low-lying en- 
ergy spectrum of the exact two-boson Hamiltonian (1) 
beyond the standard two-mode Bose-Hubbard approx- 
imation |15j . In the optical context, the scenario can 
be explained in a different view as the result of evanes- 
cent photonic tunneling among a few guided modes of 
the four two-dimensional guides in the geometrical set- 
ting of Fig. 1(b). In fact, let us indicate by 0i j2 the fun- 
damental modes of the isolated waveguides I and III, by 
03,4,5 the fundamental and the two lowest higher-order 
degenerate transverse modes of the isolated waveguide 
II, and by 06,7,8 the fundamental and the two lowest 
higher-order degenerate transverse modes of the isolated 
waveguide IV. A typical profile of such modes is shown 
in Fig. 1(c). Let then expand the envelope ip as a super- 
position of such modes with z-varying coefficients, i.e. 
tjj = J2i—xCi(z)<pi(xi,X2)exp(ij3z), where /3 is a reference 
propagation constant. Note that, for symmetry reasons, 
one has c§ — c 3 , c-j — C4 and c§ = Cg. In the tight-binding 
and nearest-neighboring approximation, neglecting cross- 
coupling terms, the following coupled-mode equations 
for the amplitudes q can be derived (see, for instance, 
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FIG. 3. (color online) (a) Behavior of pn, and (b) of P2 versus 
propagation distance, for increasing values of particle inter- 
action strength measured by the detuning <5i , as predicted by 
the coupled-mode equations (3), for «2 = 0.16, K3 = 0.80, 
and 82 = 20 (in units of mm -1 ). Curve 1: 5x = 0, «i = 0.212; 
curve 2: 61 = 1.22, ki = 0.26; curve 3: <5j = 3.2, n\ = 0.32; 
curve 4: 61 = 18.9, m = 0.38. 

El): 

i(dc\jdz) = — 2kic 3 — 2K2C4 — 2K3C5 + 5\Ci 
iidcz/dz) = —2kic 3 — 2K3C4 — 2k 2 c 5 + <5ic 2 
i(dc 3 /dz) = — «i(ci + c 2 ) (5) 
i(dc±/dz) = — k 2 Ci - K3C2 + £ 2 c 4 
i(dc 5 /dz) = -K3C1 - k 2 c 2 + S 2 c 5 

where K\ , k 2 and n 3 are the coupling constants between 
the couples of modes {0i,03}, {01,04} and {01,0s}, re- 
spectively [see Fig. 1(c)], Si — j3 3 — fti is the mismatch 
between the propagation constants /3 3 and Pi of modes 

03 and 0i, and S2 = fi 3 — Pa is the mismatch between 
the propagation constants /3 3 and (3$ of modes 03 and 

04 (or (^5). Initial condition for Eqs.(5) is q(0) = 5;,i. 
In terms of the amplitudes c/ , the percentage of bosons 
in the right well and the same-site boson probability, as 
defined by Eqs.(3) and (4), take take the simple form 

M (z) = | C i| 2 + | C 3| 2 + |c4| 2 + |c 5 | 2 (6) 

P2 M = | C i| 2 + M 2 , (7) 

respectively. In the absence of interaction, i.e. for 
An 2 = 0, one has Si — whereas 62 is much larger than 
the coupling constants. Hence, the higher-order trans- 
verse modes of waveguides II and IV are not excited, i.e. 
one has C4 ~ C5 ~ 0, and the evolution of ci, c 2 and 
C3 can be calculated exactly, yielding pr{z) = cos 2 (kiz) 
and P2(z) — (1/2) [1 + cos 2 (2ki.z)]: this is precisely the 
dynamical behavior of uncorrelated bosons (curve 1 of 
Fig. 2). As the interaction An 2 is increased, the detuning 
Si increases, whereas 82 does not change. The coupling 
constants k%, k 2 and n 3 are given by overlapping integrals 
involving the coupled guided modes, and are expected to 
slightly increase as An 2 is increased because of the less 
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FIG. 4. (color online) Light tunneling dynamics in a four- 
core fiber, (a) Profile of the optical potential n a — n(x\,X2) 
(core diameter 2w — 5 /im, core spacing 2a — 7 /im, index 
change An = 0.005). (b) and (c) show the evolution of pr 
and P2 , respectively, for increasing values of the width w c of 
the cut in guides I and III. Curve 1: w c = (non-interacting 
bosons); curve 2: w c = 0.6 /im (correlated pair tunneling); 
curve 3: w c — 2 fim (tunneling of a fragmented pair). 

confinement of modes <j>i and <j>2- If (#2 — is still large 
enough that the higher-order transverse modes of waveg- 
uides II and IV are still out of resonance, the amplitudes 
C4 and C5 remain small, and the tunneling dynamics is 
mainly governed by the first three equations of the system 
(3), but with Ji^O. A nonvanishing value of the detun- 
ing 8\ is responsible for the doubly-periodicity of Pr(z), 
the increase of the tunneling period, and the appearance 
of correlated pair tunneling (i.e. p%{z) — 1) as observed 
in curves 2 and 3 of Fig. 2. As the interaction Ari2 is fur- 
ther increased, excitation of the higher-order transverse 
modes of waveguides II and IV can not be anymore ne- 
glected, and the tunneling dynamics requires to account 
for the full five amplitudes entering in Eq.(5). For very 
strong interactions, corresponding to the fermionization 
limit, the fundamental modes 4>i t 2 of waveguides I and III 
get close to resonance with the (degenerate) transverse 
modes (j)^ an d <^7,8 of waveguides II and IV, whereas 
their fundamental modes </>3 ; 6 are now out of resonance. 
Hence, in the fermionization limit one can set c 3 ~ in 
Eqs.(5). Such equations well describe the restoration of 
tunneling of a fragmented pair. 

A typical dynamical evolution of pr(z) and p%{z) in 
the various parameter regions, as obtained by numeri- 
cally solving the coupled-mode equations (5) by varying 
<5i and taking into account for the correction of n\ solely, 
is shown in Fig. 3. Note that the behavior of both the 



percentage of bosons in the right well and the same-site 
boson probability reproduces very well the different tun- 
neling regimes previously found in Fig. 2. 

The good description of the tunneling dynamics of- 
fered by the coupled-mode equations (5) indicates that 
the tunneling dynamics of two bosons, shown in Fig. 2, is 
rather insensitive to the specific shapes of the guides, and 
could be thus observed in simpler optical structures. For 
example, in Fig. 4 it is shown that a similar dynamical 
behavior can be realized using a microstructured optical 
fiber with four circular cores of radius w and step-index 
An, in which a cut with variable width w c is applied to 
the cores I and III to mimic boson repulsion [Fig. 4(a)]. 
As the cut width w c (i.e. the interaction strength) is in- 
creased, a transition from independent Rabi oscillations 
(curve 1) to correlated pair tunneling (curve 2) and to 
tunneling of a fragmented pair (curve 3) ic clearly ob- 
served. 



IV. CONCLUSIONS 

In conclusion, an optical realization of the tunnel- 
ing dynamics of two interacting bosons in a double-well 
potential, based on light transport in a four-core mi- 
crostructured fiber, has been proposed. The present re- 
sults indicate that photonic systems could provide an 
experimentally accessible test bench to investigate in a 
purely classical setting the dynamical aspects embodied 
in the physics of strongly-correlated few-particle quan- 
tum systems. As compared to quantum simulators based 
on the coherent dynamics of cold atoms or ions trapped in 
optical lattices [TU] , the use of a classical optics simulator 
enables a direct access to the evolution of the multiparti- 
cle probability density and could provide a new route to 
realize other many-body physical models [23, 24J. For ex- 
ample, the introduction of gain and loss regions in the op- 
tical structure could offer the possibility to test in the lab 
the physics of many-body particles within non-Hcrmitian 
■pT-symmetric models j2"4"] . 
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